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We study theoretically the eftect of interactions between quantum phase slips in a short super- 
conducting wire beyond the dilute phase slip approximation. In contrast to the smooth transition 
in dissipative Josephson junctions, our analysis shows that treating these interactions in a self con- 
sistent manner leads to a very sharp transition with a critical resistance of Rq = /i/(4e^). The 
addition of the quasi-particles resistance at finite temperature leads to a quantitative agreement 
with recent experiments on short MoGe nanowires. Our treatment is complementary to the the- 
ory of the thermal activation of phase slips, which is only valid for temperature at the vicinity of 
the mean field metal to superconductor transition. This self consistent treatment should also be 
applicable to other physical systems that can be mapped onto similar sine-Gordon models, in the 
intermediate-coupling regime. 
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One of the most intriguing problems in low- 
dimensional superconductivity is the understanding of 
the mechanism that drives the superconductor-insulator 
transition (SIT). Experiments conducted on quasi-one- 
dimensional (ID) systems have shown that varying the 
resistances and diameters of thin metallic wires can sup- 
press superconductivity [l|, 0] , and in certain cases lead 
to an insulating-like behavior [1, S S] • We focus on 
recent experiments conducted on short ultrathin MoGe 
nanowires Q that show a SIT driven by the wires' normal 
state resistance with a critical resistance Rc ~ Rq- 

The universal critical resistance may suggest that at 
a temperature much lower than the mean field transi- 
tion temperature, T <ti Tc, the wire acts as a supercon- 
ducting (SC) weak link resembling a Josephson junction 
(JJ) connecting two superconducting leads. Schmid l7| 
and Chakravarty [1] showed that a JJ undergoes a SIT 
as a function of the junction's shunt resistance, Rs, 
with a critical resistance of Rq, and that the temper- 
ature dependence of the resistance obeys the power law 

R{T) oc T y ) . The theory was later extended to 
JJ arrays and SC wires B E [ill- Within these theo- 
ries, a similar power-law temperature dependence of the 
resistance at low temperature prevails. 

A power-law R{T) implies a rather wide SIT, with a 
resistance that has only a weak temperature dependence 
in the critical region. Contrary to the theoretical analogy 
with a single JJ, A. T. Bollinger et al. |6i] observe that 
the resistance of quasi ID MoGe nanowires exhibits a 
strong temperature dependence, even close to the SIT. 
In fact, the resistance could be fitted with a modified 
LAMH theory 12, l3| of thermally activated phase slips 
down to very low temperatures. Nevertheless, for these 
narrow wires it appears that the LAMH theory is valid 
only in a narrow temperature window [3,[lH|- Moreover, 
the LAMH analysis does not explain the appearance of 



the universal critical resistance Rq. 

In this manuscript we propose an approach that cap- 
tures both the universal critical resistance at the SIT, 
and the sharp decay of the resistance as a function of 
temperature. As in previous works, we treat the SIT in 
nanowires as a transition governed by quantum phase-slip 
(QPS) proliferation. This picture alone, however, cannot 
account for the observed strong temperature dependence 
of the resistance. We claim that the key ingredient left 
out in previous works is the inclusion of interactions be- 
tween phase slips in the finite-size wire, when the phase- 
slip population is dense - a situation that occurs in thin 
wires at high temperatures (see also Ref. ^llj). 

We treat these interactions in a mean-field type ap- 
proximation: when analyzing the behavior of a small 
segment of the wire, we include in its effective shunt- 
resistance the resistance due to phase slips elsewhere in 
the wire. This scheme is motivated by numerical analysis 
of a related problem, an interacting pair of resistively- 
shunted JJs [3]. Primarily, this self-consistent treat- 
ment produces the sharp temperature dependence of the 
resistance. In addition, we include the effects of the Bo- 
goliubov quasi-particles, which couple to the potential 
gradient created by each phase-slip [l3|- Consequently, 
the resistance obtained in the experiment can be fitted 
without resorting to the LAMH theory beyond its limit 
of validity. This self-consistent approximation of phase- 
slip interactions should be applicable to similar multiple 
sine-Gordon models in the intermediate-coupling regime. 

The microscopic action for a SC wire can be obtained 
from the BCS Hamiltonian by a Hubbard-Stratanovich 
transformation followed by an expansion around the 
saddle point 18|. In the limit of low energy scales. 
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u,Dq^ < Aq, this yields [19|: 

l/T 



S = N^A^l dx drf^^ [In (p^) - l] 
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where L and A are the wire's length and cross section, 
respectively, = ttD/SAo, Vp = yj {Sn/2)DAo the am- 
plitude velocity, — y^7il5Ao(2Al4^V(rTT) the phase 
velocity, Vc the Fourier transform of the short range 
Coulomb interaction, A'o the density of states, D the 
electronic diffusive constant in the normal state, and the 
SC order parameter is parameterized as A = Aope*'^, 
with Aq, the mean field solution. For the wires in Ref. |6| 
2AVcNo oc N± 1, leading to Wp <C cx Vpy/N±. Here 
the number of ID channels in the wire. 

This action supports QPS excitations, which are char- 
acterized by two distinct length scales: Vp/Ag (x ^ <C 
£,\/N± oc v^/Aq. For very long wires, ^ ^ £.VN± <C L, 
in the dilute phase slip approximation, this problem can 
be mapped onto the perturbative limit of the 1-1-1 sine- 
Gordon model In the opposite limit of very short 
wires, L < ^ <C ^VN±, the system resembles a JJ and 
can be mapped onto the 4- 1 sine-Gordon model. How- 
ever the MoGe nanowires in Ref. appear to be in the 
intermediate regime, ^ <C i <C Hence, while 

phase slips occur in different sections of the wire, they 
are indistinguishable, as each creates a phase fluctuation 
that spreads over distances larger than the wire itself. 

Moreover, the wires in Ref. [6|] have a sizable bare fu- 
gacity. Using Eq. ^ , one can estimate the core action of 
a phase slip of duration tq = I/Aq by minimizing the ac- 
tion with respect to the diameter of a phase slip . 
This yields a core action of Sc = tt/S 1/1/8(1 — e)i?Q / i?^ . 
Here e = 1 — T/Tc is the reduced temperature and i?^ — 
Rw£,/L is the normal resistance of a section of length ^. 
With this expression, the measured values of i?w and L 
(Table m, and the assumption that ^ ~ 20nm, the phase 
slip fugacity is estimated by C = e"'^'^ w 0.05 — 0.45, 
for the different wires. Namely, in the critical region, 
i?w ~ Rq, there is a dense population of phase slips 
which interact with one another; thus the dilute phase 
slip approximation is no longer a proper description. 

The flow equation for the fugacity of a phase slip any- 
where in the wire is given by: 
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where dl — —din A, and A is the running RG scale. Eq. 
([2]) treats the phase slip as occuring on an effective J J, 
with Rs being the effective shunting resistance of the en- 
tire wire. When ( is small, Rs will include only the ef- 
fective impedance of the leads. If ( is not very small, we 
will need to include in Eq. ^ additional terms of higher 
powers of C, which describe interactions between QPSs. 



Slips, Hps - and the quasi- 
[Hl, Rau- Once the distur- 



Currently, there is no complete understanding of how to 
approach the field theory of phase slips in the regime of 
intermediate fugacity C ^ 1- Therefore, we suggest to 
take into account a finite ^ by including the resistance 
due to other phase slips in the wire in the effective shunt 
resistance of the junction, R^. This treatment was suc- 
cessfully tested by numerical analysis of a simpler analog 
of the system, an interacting pair of resistively shunted 
JJs [3]. Indeed, including the interaction between phase 
slips as a modification of the shunt, Rs, is akin to guess- 
ing the form of a complete resummation of higher-order 
C terms in Eq. ([2]) [H, HI . 

The main physical intricacy is the determination of the 
effective shunt resistance, i?s(C)j that governs the rcnor- 
malization of the phase slip fugacity [Eq. ^] . A phase 
slip produces time- varying phase gradients, and hence 
electrical fields. These dissipate through two channels 
in parallel: the SC channel - which has an effective re- 
sistance due to other phase slips, Ru 
particles conducting channel ij^jj, ±i.qp. 
bance reaches the leads, it also dissipates through the 
electro-dynamical modes of the large electrodes, whose 
real impedance is parameterized by i?cioc- 

For T ^ Tc, the resistance of the quasi particles, Rqp, 
can be approximated by Rqp = e^r"nc,p i " ~ 

Rn\ 27rAo ^~ ' Unfortunately, we lack a microscopic 
model for the impedance of the electrodes, as this de- 
pends on the details of the system. However, we expect 
that at large scales, T < A < A, the electrodes will act 
as a transmission line to the electromagnetic waves gen- 
erated by the phase slip. This transmission line is char- 
acterized by a real impedance which we denote as i?eiec, 
and use as a fitting parameter. Hence, the effective shunt 
resistance that affects the renormalization of the phase 
slip fugacity at T < A < A [Eq. (O] is 

i?s[C(A)] = i?elcc + 1/ (l/i?ps[C(A)] + l/RqpiT)) . (3) 

Fig. [1] shows the circuit we suggest describes the system. 

The resistance is measured in response to an applied 
DC current. At this zero frequency limit, the electrodes 
act as a capacitor connected in parallel to the wire. 
Therefore, the measured resistance is the total wire re- 
sistance, unaffected by the environment, which is cut off 
from the wire. 

i?tot(T) = 1/ (l/i?ps[C(T)] + l/i?qp(T)) . (4) 

The occurrence of a phase slip causes a resistance in the 
otherwise SC wire through the relation Rps oc (L/^)C^. 
Using this relation and Eq. ([2]), we can write an RG 
equation for the dimensionless resistance 

^ = 2(^.2^1-^ 

dl ^ dl V -Rs(C) 

diRps/Rq) ^2(1- — ^ [ 

dl \ Rs{Rps)J \Rq 



(5) 
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FIG. 1: The effective electrical circuit for the nanowire. A 
phase slip produces time-varying phase gradients, and hence 
electrical fields, that dissipate through the quasi-particles con- 
ducting channel, -Rqp, and through the SC channel, that has 
an effective resistance due to other phase slips, Rps. Once the 
disturbance reaches the edges, it also dissipates through the 
electro-dynamical modes of the large electrodes, whose real 
impedance is represented by Rdcc- 



If we now substitute the effective resistance in Eq. ([3]), 
depicted in Fig. [1] into the RG equation (O, we find the 
following self consistent equation: 



dl 



= 2 



1 - 
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Integrating Eq.® from the ultra violet (UV) cutoff, 
A(T*) = T*, to the infra red cutoff, T, yields i?ps(r). 

The wire's DC resistance Eq. (U)), calculated using Eq. 
([6]) , is plotted in Fig [2] as a function of temperature for 
different R*. Here R* = Rtot{T*) is the normal state 
resistance of the wire, at the UV cutoff A(r*) = T*. 
We assume that the wires are thin enough such that the 
mean field transition from normal to SC is wide, and 
R* « RtotiT^) = Rn- For simplicity, we have assumed 
throughout our calculation that A(r) w A. This as- 
sumption holds at the low temperature regime, T <C A, 
where the theory of QPSs, based on the effective action 
Eq. ID), is expected to be vahd. The resistance of the 
environment is taken to be i?eiec = O.IRq. In practice, 
A, R* and i?cioc can be used as fitting parameters, when 
comparing the theory to the experiment. Moreover, R* 
can be determined independently as the resistance mea- 
sured below the drop that indicates passing through Tc 
of the SC films (see Ref. S). 

Eq. ([6]) is also applicable to a wire with a normal re- 
sistance larger than the quantum resistance, R* > Rq. 
In this limit the fugacity increases in the renormaliza- 
tion process, and Eq. ([2]) is no longer valid for ( > 1. 
For wires with R* > Rq, we overestimate the phase 
slip fugacity as ({T*) — 0.5, and plot the reduced resis- 
tance, Rtot/Rg, down to the temperature, Tq, for which 
C(To) — 1. The results are shown in Fig. [2] Overestimat- 
ing C{T*) gives an upper bound on Tq, where Eq. ^ 
is no longer applicable. Fig. [5] shows that the transi- 
tion between SC and insulating wires occurs for a criti- 
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FIG. 2: Reduced resistance Rtot/Rg as a function of renor- 
malized temperature T/A for different R* /Rq, ranging from 
0.8 (lower plot) to 1.01 (upper plot) in increasing steps of 
0.03. The environment resistance was taken to be Reiec/Rq = 
0.1. Gray traces indicate the total wire resistance _Rtot = 



-f?pB-f?qp(r) 



with Rps{T) = R* (T/A) 



2{Rq/R,-1) 



and con- 



flps+-Rqp(T) ' 

slant Rs = K* -\- Rdcc- The insulating plots are cutoff at 
temperature. To, for which C,{To) = 1. 



cal resistance Rc ~ Rq. But in contrast to the standard 
Josephson junction theory, (gray curves in Fig. [2]), the 
transition is much more sharp (notice the log scale). 
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177 


5.46 


4.2 


2.5 


1.2 
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43 


3.62 


2.62 


2.35 


1.25 
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63 


2.78 


2.13 


3.07 


0.66 
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93 


3.59 


2.89 


3.85 


0.55 


5 


187 


4.29 


4.5 


6.55 


0.31 


6 


99 


2.39 


2.09 


4.84 


0.4 



TABLE I: Summary of nanowire parameters, and the param- 
eters used to fit the experimental data; L is the length of the 
wire determined from SEM images, Rw is the wire's normal 
state resistance, taken as the resistance measured below the 
film transition. Fitting parameters: R* = Rtot {T — A) is the 
wire's resistance at the UV cutoff, A the SC order parameter, 
and -Rcicc the impedance of the electrodes at T < A < A. 

A comparison between the theoretical curves and the 
experimental data taken from Ref. is shown in Fig. [31 
The curves were calculated by fitting R*, A, and i?cicc- 
Since the theory of QPSs is expected to be valid at 
T <C A, deviations from the theoretical curves at high 
temperature are reasonable. In general, as A is propor- 
tional to Tc, increasing A shifts the sharp decay of the 
resistance to high temperatures. Both R* and i?eioc af- 
fect the high temperature resistance, whereas R* and A 
control the width of the transition. 
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We have made an attempt to fit the experimental data 
corresponding to the insulating wires of Ref. |6|. As the 
insulating wires are thinner, we expect a strong suppres- 
sion of Tc 0- Consequently, the theory of QPSs with 
^ < 1, is valid in a narrow range of temperatures. While 
we manage to fit the experimental curves in this regime of 
parameters, we do not find it of scientific merit to present 
the results (our fits cover a range of ~ 10 data points, 
with three fitting parameters). 
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FIG. 3: Comparison between the theory (solid line) and the 
experimental data [^. Details of the wires, as well as fitting 
parameters are summarized in Table HI Note that the esti- 
mated values for R* are in good agreement with the measured 
values of the resistance of the wire at T = A. The discrepan- 
cies between theory and experiment occur for the longer wires 
(curves 1 and 5). In these wires the condition v^/Aq 3> L 
might not hold, and one should consider additional renormal- 
ization of the fugacity at higher scales. Inset: The effect of 
an adjustment of the fit parameters, in roughly 10% (from 
dark gray to light gray): A = 4.356Jf, Reiec = 4400 and 
R* = 188m. 



In conclusion, we have studied theoretically the effect 
of interactions between QPSs, in short SC wires, beyond 
the dilute phase slip approximation. Our analysis shows 
that treating these interactions in a self consistent man- 
ner, produces a sharp SIT transition, with a critical re- 
sistance Rc ~ Rq, in agreement with recent experiments 
conducted on short ultrathin MoGe nanowires i6| . More- 
over, we have shown that adding the resistance due to 
the presence of a finite number of BdG quasi-particles 
in the wire leads to a quantitative agreement with the 
experimental curves. Our method should be applicable 
to a wider range of physical problems which involve the 
proliferation of topological defects with a sizable bare fu- 
gacity. In particular, it could be applied to the study of 
a Luttinger liquid with an extended impurity [l^l . 
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